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Semileptonic Decays of D and B Mesons* wub 95-34 
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We report results of our ongoing investigation concerning semileptonic decays of heavy pseudoscalar mesons 
into pseudoscalar and vector mesons. Particular attention is paid to uncertainties in the q 2 and the heavy quark 
mass dependence of formfactors. Moreover we present a non-perturbative test to the LMK current renormalization 
scheme for vector current transition matrix elements and find remarkable agreement. 
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1. INTRODUCTION 

The accurate determination of Kobayashi 
Maskawa (KM) matrix elements involving heavy 
quarks requires control over the low energy QCD 
parts of the corresponding weak transition ampli- 
tudes. Although lattice QCD is the per se method 
to tackle this regime, it suffers from fact that the 
lattice resolution is still too poor to represent b 
quarks directly on the lattice. 

With our current high statistics project we 
push forward to a lattice resolution of a -1 ~ 3.2 
GeV. This allows for direct calculation of transi- 
tion amplitudes up to a mass region of (1 — 1.6) x 
m£>, keeping discretization errors small. As direct 
access to the b mass is still excluded, we work with 
four heavy quark masses, which provide us with 
a sufficient lever arm for the extrapolation to the 
B meson. 

2. LATTICE SETUP 

We are working on 24 3 x 64 lattice at (3 = 6.3, 
with standard Wilson quarks in the quenched ap- 
proximation. The heavy quark masses are rep- 
resented by K h = 0.1200,0.1300,0.1350,0.1400. 
This covers the physical region 0.8m c < 
rrih < 1.6m c . For the light quark masses we 
have chosen kj = 0.1450,0.1490,0.1507,0.1511, 
which corresponds to 0.8m 3 < m; < 3m,. 
The decay amplitudes have been constructed 
with the initial meson at rest and the fi- 



nal meson carrying spatial momenta pa = 
|| x {(0,0, 0); (1,0,0); (1,1,0); (1,1,1); (2, 0,0)} 
+ permutations. In order to improve on the 
groundstate projection we have used Wuppertal 
"gaussian" smeared wavefunctions for the quark 
fields. Details can be found in ref.[?]. 

The simulation has been performed on the 32 
node CM5 located at the IAI in Wuppertal. We 
have created a total of 100 independent gauge 
configurations, which served for calculation of 2- 
and 3-point correlators. 

The analysis has been done so far on a sub- 
set of 60 configurations. Using this data set we 
find n„it = 0.151818(37). The lattice spacing 
has been determined with three different meth- 
ods. With R. Sommer's [?] method we find 
= 3.314(29) GeV, extrapolation of the vec- 
tor meson mass to n cr i t yields o" 1 = 3.453(217) 
GeV, and with the stringtension ^fa = 440 MeV 
we get a" 1 = 3.210(52) GeV. The results agree 
nicely within errors. In the following we will use 
URo to set the scale. 



3. q 2 

The KM matrix elements can be extracted from 
experimental measurements of the decay widths 
r and the branching ratios R. In doing so one 
needs the formfactors 2 F(q 2 ) , which parametrize 
the transition matrix elements[?]. As the major 
contribution to the phase space integrals for T 
and R comes from the region of small q 2 , one con- 
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2 For the time being F(q 2 ) denotes generically the form- 
factors /+ , f , V, Ai , A 2 . 
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veniently quotes the values F(q 2 = 0). Of course 
knowledge of F(q 2 ) over the entire q 2 region is 
neccessary in order to calculate R and T accu- 
rately. 

According to our setup of spatial momenta, 
the calculation provides us with (up to) 5 entries 
F(q 2 ) for each mass combination, which we use 
to extrapolate to -^(O) as well as to learn about 
the functional dependence of F(q 2 ). 

In our analysis we have tested three different 
methods to fit our data: Two of them are guided 
by the pole dominance hypothesis 

nf) = ^r, (1) 

where we a) took a fixed pole mass m ( , as ex- 
tracted from the corresponding 2-point correla- 
tors, and b) treated m ( as a free parameter. The 
third Ansatz assumes linear behaviour 

F(q 2 ) = a + bq 2 . (2) 

Fig. 11 shows a typical example of our formfactor 
data for weak transitions of a heavy-light pseu- 
doscalar meson into light-light pseudoscalar and 
vector mesons (HI — » 11), together with the dif- 
ferent fits. It turns out that all Ansatze describe 
the data reasonably well and lead to compatible 
results at q 2 = 0. Therefore the good news is that 
we are able to determine the value of F(q 2 = 0) to 
come out independent of the formfactor model. 

On the other hand, however, our data does 
not clearly distinguish between the different func- 
tional forms. This is due to the statistical noise of 
the data and gives rise to systematic uncertain- 
ties. 

A similar situation is found for the decays of a 
heavy-light pseudoscalar meson into heavy-light 
pseudoscalar and vector mesons (HI — » H'l). 
Again, the results at q 2 = are almost unal- 
tered by the fit method, but sizeable uncertainties 
concerning the functional form of F(q 2 ) over the 
entire q 2 range remain. 

4. RENORMALIZATION 

In order to convert the above lattice data into 
continuum results, we need the renormalization 




Figure 1. q 2 dependence of formfactor s for HI 
— » 11 transitions. The mass of the heavy initial 
quark corresponds to K{ = 0.1350 and the masses 
of light final and the light spectator quark is given 
by ki = K spec = 0.1490. The solid line represents 
the fit to eq.11 with fixed m t , the dashed line keeps 
m t free and the dotted line correponds to a fit to 
eq.11. 



constants Zy and Za of the vector and axialvec- 
tor currents. It has been pointed out by C.W. 
Bernard [?] that the renormalization constant of 
the local vector current Zy depends strongly on 
the mass of the quarks involved. Following his ar- 
guments this undesired mass dependence can be 
removed if the (standard) \/2k normalization of 
the quark fields is replaced by the LMK[?] pre- 
scription. 

As a first step to check this issue 3 with our data 
we study the ratio 

R n _ (P(m^,mQ|y — \P(m H , mi )) f 
H ' H{P > {P{ m ' H , mi )\V^\P{m H , mi )) f - ( > 

yioc(con$) . g ^ e o'th component of the lo- 
cal(conserved) vector current, and P denotes a 
pseudoscalar eigenstate. 

Rhh(Q) yields directly the renormalization fac- 
tor of Vq oc , as Vq 071 ' is unaltered in that instance. 
For the study of decays we are rather interested 

3 see also [?]. 
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in the ratios -R/jh(0) and Rhh(P) however. The 
LMK predictions for these ratios read 



Rhh(0) 



R hH (0) 
R hH {v) 



3 kh 



2k h 



Z l j c 1 



3 kh 



peons 
RhH(0) ■ 



2k h 



(4) 

(5) 
(6) 



Taking the standard normalization one gets 
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Figure 2. The ratios Rhh(P) (left) and Rhh(Q) as 
a function of the initial (k{) and final (k} ) quark- 
mass. The solid line represents the LMK predic- 
tion, the dashed line corresponds to the standard 
scheme. 



meson mass. As the exact mass dependence of the 
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Figure 3. Mass dependence of form] actors at 
qr 2 = . Crosses denote the data. The solid 
line refers to fit method a), the dashed line to b) 
and the dotted line to c). The extrapolated results 
are depicted at 1/ttxb as filled circles (method a), 
upper triangles (method b), and lower triangles 
(method c). 



RhH = Rhh = Z l ™. In fig. ?? we display 
the measured ratios Rhh(P) and Rhh(Q) together 
with the LMK and standard predictions. We 
have inserted Z l j c = 1 - ^-g 2 , Z l j c = 1 - 
0.174081c/ 2 , and Z^ ons = 1, using the boosted 
coupling g-\l ) = g^il/STvP^) + 0.02461. 
The LMK prediction describes the data remark- 
ably well within the entire mass range, whereas 
the standard prediction fails drastically. In the 
following we will therefore apply the LMK scheme 
to our formfactor results, although a similar ratio 
test for the spacelike components is still missing. 

5. MASS DEPENDENCE 

The determination of the formfactors for B me- 
son decays requires extrapolation in the heavy 



formfactors is still unknown we have used several 
Ansatze for the extrapolation, in order to check 
on systematic uncertainties: 

f °) a +^7 
^P5,9 2 =0)= b) 7 J=( a +^) (7) 

In fig.?? we display the formfactors F(q 2 = 0) for 
HI — >■ 11 transitions, as a function of 1/mps- Ob- 
viously the Ansatze (eq.'s ??) lead to consistent 
results at mps = vtlb- 

We mention that the transitions HI — >■ H'l ex- 
hibit a much weaker mass dependence, and again 
the results at the B mass are unaltered by the fit 
method. 
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Table 1 

Form) "actors at q 2 =0 for weak decays into pseu- 
doscalar mesons. 



mode 


/+(0) 


/o(0) 


Z) — ?■ 7T 


0.68(13)1/ 


0.65(5)1? 




0.71(12)1/ 


0.67(6)1? 


B — ?■ 7T 


0.50(14)1* 


0.20(3)1?. 


B -> D 


1.01(60) 


0.88(29) 


B S ^D S 


0.99(42) 


0.85(24) 



6. RESULTS 

In tables ?? and ?? we quote our results for 
the formfactors at q 2 = 0, which have been ob- 
tained after extra- and interpolation of the heavy 
quark mass to the b and c quark and extra- and 
interpolation of the light quark mass to m q = 
and m q = m st range respectively. The errors in- 
clude both statistical (in brackets) and system- 
atic uncertainties 4 . The results for D decays are 
in good agreement with experiment and other 
lattice calculations (see e.g. [?]), giving confi- 
dence to our analysis method. For B — > ir, p de- 
cays however, the situation is much less settled. 
The data obey the condition /+(0) = /o(0) on a 
2cr level only, indicating systematic uncertainties 
arising from the extrapolation in the heavy mass. 
A careful study of this issue will be included in 
the analysis with full statistics. 



Table 2 



Formfactors at q 2 =0 for weak decays into vector 
mesons. 



mode 


V(0) 


^(0) 


^(0) 


D -> p 


1.31(25)111 


0.59(7)1« 


0.83(20)1* 2 


D — ► K* 


1.34(24)1^ 


0.61(6)1? 


0.83(20)1* 2 


B -> p 


0.61(23)1! 


0.16(4)1^ 


0.72(35)1/ 


B -> D* 


1.84(1.10) 


0.97(32) 


2.76(1.58) 


B,^D* 


1.84(1.10) 


0.91(24) 


2.25(1.12) 
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4 This being an intermediate analysis, we have omitted in 
some cases to quote systematic uncertainties. 



